Critical exponents of domain walls in the two-dimensional Potts model 
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We address the geometrical critical behavior of the two-dimensional Q-state Potts model in terms 
of the spin clusters (i.e., connected domains where the spin takes a constant value). These clusters are 
different from the usual Fortuin-Kasteleyn clusters, and are separated by domain walls that can cross 
and branch. We develop a transfer matrix technique enabling the formulation and numerical study 
of spin clusters even when Q is not an integer. We further identify geometrically the crossing events 
which give rise to conformal correlation functions. This leads to an infinite series of fundamental 
critical exponents hi x -i 2 ,u X i valid for < Q < 4, that describe the insertion of t\ thin and ^2 thick 
domain walls. 



PACS numbers: 64.60.De 05.50+q 

Many geometrical features of two-dimensional (2D) 
critical phenomena are by now under complete control, 
thanks to the combined powers of Conformal Field The- 
ory (CFT) and Schramm-Loewner Evolution (SLE) [HE]. 
This situation is epitomized by the Q-state Potts model 
for which the Fortuin-Kasteleyn (FK) expansion of the 
lattice model gives rise to a formulation in terms of clus- 
ters and their surrounding loops (hulls). These loops 
behave like the SLE trace in the continuum limit [2] , and 
viewing them as contour lines of a (deformed) Gaussian 
free field leads to the Coulomb Gas (CG) approach to 
CFT jU [3]. Our understanding of the critical proper- 
ties of FK clusters and loops can be considered almost 
complete. 

And yet the properties of spin clusters themselves — 
i.e., the connected domains with a constant value of the 
Potts spin — have as a rule remained ill understood. This 
is particularly frustrating, since those are the very clus- 
ters that one would observe in an actual experiment on a 
magnetic alloy in the Potts universality class. The case 
of the Ising model Q — 2 is an exception to this rule 
[U [5] , but this is due to its "coincidental" equivalence to 
the 0(n) vector model with n = 1. Indeed, denning the 
Ising spins on the triangular lattice, the corresponding 
O(n) model is described by self and mutually avoiding 
loops on the hexagonal lattice, and such loops are read- 
ily treated by CFT and SLE techniques. 

For general Q, the salient feature of Potts spin clus- 
ters is that the domain walls separating different clusters 
undergo branchings and crossings (see Fig. [I]). These 
phenomena are however absent for Q = 2 with the 
above choice of lattice. It is precisely these branch- 
ings and crossings that make the application of exact 
techniques — such as CG mappings or direct Bethe ansatz 
diagonalization — very difficult, if not impossible. The be- 
lief that spin clusters are indeed conformally invariant for 
other values of Q ^ 2 in the critical regime < Q < 4 
has even been challenged at times, but seems however 




FIG. 1: A configuration of the Q — 3 Potts model, and the 
corresponding set of branching domain walls. 

well established by now [5J [7] . 

Some progress has been accomplished in the Q = 3 
case H] by speculating that the spin clusters in the 
critical Potts model would be equivalent to FK clusters 
in the tricritical Potts model |9]. This equivalence 
has however not been proven, and is moreover restricted 
so far to the simplest geometrical questions [ID] . The 
equivalence can also be understood as a relationship with 
the dilute O(n) model [II]. 

The Potts model has been used recently to build a new 
class of 2D quantum lattice models that exhibit topo- 
logical order [12]. Both FK clusters and domain walls 
between spin clusters are important in these models. 

Apart from issues of branching and crossing, another 
major hurdle in the study of Potts spin clusters has come 
from the lack of a formulation that can be conveniently 
extended to Q a real variable. In the case of FK clus- 
ters, this formulation led naturally to the introduction 
of powerful algebraic tools via the Temperley-Lieb (TL) 
algebra, and to the equivalence with the 6- vertex model — 
the eventual key to the exact solution of the problem [1] . 
Factors of Q then appear naturally through a parameter 
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in the TL algebra, or — via a geometrical construction — 
as complex vertex weights in the 6-vertex formulation. 
Also for spin clusters can Q be promoted to an arbi- 
trary variable: the weight of a set of spin clusters is 
simply the chromatic polynomial of the graph dual to 
the domain walls. From the point of view of the TL al- 
gebra, the domain walls are composite (spin-1) objects, 
hence more complicated than the FK loops. Recent work 
on the related Birman-Wenzl-Murakami (BWM) algebra 
[131 H3] suggests that this formulation might be amenable 
to the standard algebraic and Bethe ansatz techniques, 
although such a lofty goal has not been achieved so far. 

We report in this Letter major progress towards the 
understanding of Potts spin clusters. Our results are 
of two kinds. On the one hand, we develop a transfer 
matrix technique which allows the formulation and nu- 
merical study of the spin clusters for all real Q. On the 
other hand, we identify the geometrical events that give 
rise to conformal correlations, and provide exact (albeit 
numerically determined) expressions for an infinite fam- 
ily of critical exponents, similar to the familiar "L-legs" 
exponents [TJ [3] for FK loops. Surprisingly, we find that 
geometrical properties of spin clusters encompass all in- 
teger indices (r, s) in the Kac table h ryS . An analytical 
derivation of our results appears for now beyond reach, 
in part because the algebraic properties of our transfer 
matrix are still ill understood. We do however provide 
some exact results based on an approach which does not 
involve a CG mapping, but rather the use of a massless 
scattering description. 

Domain wall expansion. The Q-state Potts model is 
defined by the partition function 

Z = Y, U exp(K5 ai , aj ) , (1) 
<t {ij)eE 

where K is the coupling between spins <7j = 1,2, ... ,Q 
along the edges E of some lattice C (we use the square 
lattice in the computations below, but the triangular lat- 
tice in the figures). The Kronecker delta function 5 aurTj 
equals 1 if <r, = aj, and otherwise. 

The domain wall (DW) expansion of involves all 
possible configurations of domain walls that can be drawn 
on the dual of C (see Fig. [2]). A DW configuration is 
given by a graph G (not necessarily connected). The 
faces of G are the spin clusters. Since we do not specify 
the color of each of these clusters, a DW configuration 
has to be weighted by the chromatic polynomial Xg(Q) 
of the dual graph G. Initially Xq(Q) ^ s defined as the 
number of colorings of the vertices of the graph G, using 
colors {1, 2, . . . , Q}, with the constraint that neighboring 
vertices have different colors. This is indeed a polynomial 
in Q for any G, and so can be evaluated for any real Q 
(but Xg(Q) is integer only when Q is integer). For exam- 
ple, the chromatic polynomial of the graph G on Fig. [2] is 
Q(Q - 1) 7 {Q - 2) 7 . The partition function ([!]) can thus 




FIG. 2: A domain-wall configuration corresponding to a graph 
G, and its dual graph G. 

be written as a sum over all possible DW configurations 
Z = o NK J2(-- K ) lCagtHG) Xa(Q) (2) 

G 

where N is the number of spins, and length(G) denotes 
the total length of the domain walls. 

The bulk DW exponents. The fundamental geomet- 
ric object we consider is a connected part of a domain 
wall that separates two clusters. One can ask how the 
probability, that a certain number of such DW connect a 
small neighborhood A to another small neighborhood B, 
decays when the distance x between A and B increases. 
Each DW separates two spin clusters which connect A 
and B. There are in fact two types of such DW, de- 
pending on the relative coloring of the two clusters that 
are separated. If the two clusters have different colors, 
they can touch, so the DW is thin (see Fig. [3]a). If the 
two clusters have the same color, then they cannot touch 
(otherwise they would not be distinct), so the DW has 
to be thick (see Fig. [3jb). 

We can now state the central claim of this Letter. Con- 
sider the 2D Potts model for any real Q in the criti- 
cal regime < Q < 4. Then the probability P that 
the two regions A and B, with separation x 3> 1, are 
connected by i\ thin DW and £2 thick DW decays al- 
gebraically in the plane with some exponent h{Q,t\,l-i), 
viz., P oc a^MQAA). Equivalently, on a long cylinder 
of size L X I with I ^> L, and A and B identified with the 
opposite ends of the cylinder, the decay is exponential: 
P oc e ~ 4 ^( e / L ) h (Q^i^2) _ Below we check this assertion 
numerically, and we observe that the numerical values of 
the exponents match the formula 

h(Q,£ 1 ,t 2 ) = h^-i^n-L , (3) 

where we have used the Kac paramctrization of CFT 

hrs= (r- SK /4f-(l- K /4f 

K 

and 2 < k < 4 parametrizes Q = 4 (cos <E [0,4]. 
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surrounding a single lattice edge. This edge links spins 
(shown as solid circles) on diametrically opposite sites of 
the rhombus. On an L = 4 square lattice with periodic 



boundary conditions this reads T ■ 

A rhombus corresponds to a vertical edge, and acts 
on a basis state s as follows. If exactly Q s distinct color 
labels {ck} are used in s, then the new color label c[ of 
the spin o~i can be either unchanged c\ = (with weight 
e K ), or any one of the other labels already in use c- = c k 
(each with weight 1), or a new one c[ (jt {c k } (with weight 
Q — Q s ). Note that this latter weight is in general non- 
integer, and is responsible for the correct computation of 
the chromatic polynomial Xq(Q)- A rhombus "0* adding 
a horizontal edge between vertices i and i + 1 corresponds 
simply to a diagonal matrix, with a weight e K if Cj = Cj+i, 
and 1 otherwise. 

With these rules at hand, one can write the (periodic) 
L = 3 transfer matrix for arbitrary Q in the basis ^ as 
an instructive example: T — hi • hi • /13 • i>i • i>2 ■ v 3 with 



FIG. 3: The two different types of DW: a thin DW corre- 
sponds to the interface between two clusters of different col- 
ors (a), while for a thick DW the two clusters have the same 
color. An illustration for the Q — 3 Potts model is given (left) 
as well as a schematic picture for non-integer Q (right). 
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Transfer matrix formulation. The DW expansion ^ 
may appear unwieldy and difficult to study in a Monte- 
Carlo simulation for non-integer 0- However, it can be 
tackled in a transfer matrix formalism that is no more 
complicated than the one |15j routinely used in the study 
of the FK clusters. 

Consider a strip of the square lattice of width L spins 
(boundary conditions will be detailed later). The basis 
states on which the transfer matrix T acts contain one 
color label Ci per spin. By definition, one has Ci = Cj 
if and only if <7j = o~j (i.e., the two spins on sites i 
and j have the same color). The color labels c; contain 
less information than the spin colors o~i themselves. For 
instance, any configuration in which the first and third 
spins have the same color, no matter which one, and no 
other spins have identical colors, is represented by 



Ci c 2 Ci C4 Cl 

• • • • • • • • 



(5) 



For a row of L — 3 vertices, and for any (non-integer) Q, 
there are thus precisely five basis states 



Ci Ci Ci 

• • • 



Cl C 2 Ci 

• • • 



Cl Ci c 2 
• • • 



Cl c 2 c 2 
• • • 



Cl c 2 c 3 

• • • (6) 



Note that the last state would carry zero weight for Q = 
2, but apart from that the number of basis states for any 
given L will be finite and independent of Q. 

We can write T as a product of elementary transfer 
matrices, each represented symbolically as a rhombus 



and h\ = diag(e , l,e , 1, 1). The remaining matrices 
can be obtained from those given by cyclic permutations 
of the color labels. The reader can now check that T 
gives the same free energy as the FK transfer matrix, 
even when Q is non-integer. 

However, to obtain the desired two-point correlation 
functions of DW, the basis states Q need to be endowed 
with some additional information about the connectivity 
of the spin clusters. We need to know whether two spins 
having the same label Ci also belong to the same cluster. 
Thus the states we use in the final transfer matrix have 
the form 



Cl c 2 Ci Cl c 5 
• • • • • 



V 



ci c 2 ci ci c 5 



(T) 



In the left state, the spins on vertices 3 and 4 are in the 
same cluster, but not in the same cluster as the spin 1. 
In the right state the spins 1, 3 and 4 are all in the same 
cluster. These two states are different. In the transfer 
matrix evolution, each time two neighbor vertices corre- 
spond to the same color, the corresponding clusters are 
joined up. 

The final transfer matrix thus keeps enough informa- 
tion, both about the mutual coloring of the sites and 
about the connectivity of the clusters, to give the cor- 
rect Boltzmann weights to the different configurations, 
even for non-integer Q, and to follow the evolution of the 
boundary of a particular set of clusters. These bound- 
aries are precisely the domain walls (see Fig. [3]) . 
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2 
3 
4 
5 

Exact 



(2,0) (0,2) (3,0) (2,1) (0,3) 



2.01(1) 5.99(2) 2.97(4) 8.94(2) 

4.01(1) 7.99(2) 6.02(2) 8.04(2) 11.98(3) 

5.93(2) 10.01(2) 8.89(5) 10.93(5) 15.05(4) 

7.77(4) 12.09(8) 11.6 (1) 13.8 (1) 18.2 (2) 

2(p - 1) 2(p + 1) 3(p - 1) 3p - 1 3(p + 1) 



TABLE I: Critical exponents corresponding to five different 
DW configurations (£1,12), as functions of the parameter p = 
jz^:, along with the conjectured exact expression 030 • The 
table entries correspond to \p\, when |4| is rewritten as h r ,s = 
(p 2 — l)/(4p(p + 1)), with error bars in parentheses. 



Numerical results. We have numerically diagonalized 
the transfer matrix in the DW representation for peri- 
odic strips of width up to L = 11 spins. We verified 
that the leading eigenvalue in the ground state sector co- 
incides with that of the FK transfer matrix, including 
for non-integer Q. As to the excitations, we explored 
systematically all possible coloring combinations for up 
to 4 marked spin clusters, for a variety of values of the 
parameter k. Finite-size approximations of the critical 
exponents h were extracted from the leading eigenvalue 
in each sector, using standard CFT results pp. Final re- 
sults (see Table [l| for the exponents were obtained by 
extrapolating those approximants to the L —¥ 00 limit. 
Note that ^ is only valid for (£i,£ 2 ) = (0, 1) if the spin 
cluster is forbidden to wrap around the periodic direc- 
tion. Without that restriction we obtain /io,i/2- 

Analytical results. While it is far from obvious to de- 
rive ^ by CG methods, minor progress can be achieved 
using a rather different set of ideas. Indeed, we can learn 
about the dynamics of DW in the critical theory by using 
known information about the low-temperature (K > K c ) 
phase of the Potts model. Albeit non-integrable on the 
lattice, the corresponding deformation by the operator 
$21 is integrable in the continuum 16J. It can be de- 
scribed using a basic set of kinks K a t, separating two 
vacua, i.e., ordered regions where the dominant value of 
the spin is a, resp. b. These kinks scatter with a known 
S-matrix related to the BWM algebra [14] . Importantly, 
the dynamics conserves the number of kinks: the process 
KabKbc K ac is forbidden (as in any elastic relativistic 
scattering theory), although kinks do appear as bound 
states in kink-kink processes. Many properties of these 
kinks can be calculated using integrability techniques. 
When the mass m (i.e., K —> K c ), the S-matrix 
provides a "massless scattering" description [T7] of some 
of the degrees of freedom of the critical theory itself. It 
is not entirely clear what a kink, which is well-defined 
for K » K c , becomes at K c , but it is natural to expect 
that thick DW are described by the propagation of two 
(or more) kinks such as K a bKb a . As for thin DW, the 
potential existence of regions where they are reduced to 
a single edge — that is a single kink — suggest they have to 



do with more complicated processes involving two kinks 
merging into one. We will not discuss them here. It 
is an easy exercise to obtain the scaling dimension of 
thick DW using the massless scattering description. In- 
deed, the fact that the S-matrix satisfies relations from 
the BWM algebra allows us to reexpress it in terms of 
the a 2 or Bullough-Dodd S-matrix 18 , for which the 
thermodynamic Bethe ansatz was studied in 19J. This 
describes as well the dynamics of the field theory 
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S = 



8tt 



d 2 x 



with /3 2 = ~. Giving each kink the fugacity Q — 1 pro- 
duces the correct central charge c = 1 — | ^~ 4 ^ , where 

each kink has a U(l) charge equal to 0, ±1. The scaling 

■2 

of the sector with charge j produces a gap Aj = 4-, so 
the leading dimension is Aj — (1 — c)/24 = hj/ 20 . This 
agrees with ft._£ 2 ,o with 1% — j/2, so there are two kinks 
per thick DW. 

Fractal dimensions. The dimension of a spin clus- 



(8+k)(8+3k) 
32k 



ter is d = max(2, 2 — 2/10,1/2) = max ^2, 
in agreement with [?J [5]. The boundary of a cluster 
has dimension c?b = 2 — 2h-ifl = 1 4- |. According 
to j3j|, the set of points where a thin (resp. thick) DW 
has minimal width (one resp. two lattice spacings) has 
dimension d x = 2 — 2h 2A = (4 — k)(5k — 4) resp. 
d 2 = min(0, 2 - 2/i_ 2 ,o) = 0. The fact that d x > and 
e?2 = validates the epithets 'thin' and 'thick' for the 
scaling limit. For Q = 4 (or k = 4) we have d\= d 2 = 0, 
so thin and thick DWs become indistinguishable. Indeed 
one has hg 1 _g 2 2e 1 — {l\ + ^) 2 /4 in that case. For the 
Ising model Q — 2 (or k — 3) the absense of branchings 
means that one thick DW equals two thin DW. Indeed 
h^-es&i = (4^ 2 - l)/48 with I = l x + 2l 2 in that case, 
and this agrees with the exponent f° r ^ loop strands 

[3 [3] in the dilute 0(1) model. 
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